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Abstract 

The  movement  of  wafer  and  Bf  was  measured  in  unsaturated  and  partially 
frozen  soil  columns  subjected  to  linear  temperature  fields.  Both  water  and 
Bf  moved  from  fhe  warmer  fo  cooler  parts  in  the  columns.  The  data  were 
analyzed  under  the  assumption  that  Bf  is  confined  to  unfrozen  water.  It  was 
found  fhatthe  negafive  adsorption  of  Bf  by  clay  surfaces  plays  a  significant 
role  in  the  transport  of  Bf  in  frozen  Morin  clay  and  that  Bf  tends  to  move 
faster  than  unfrozen  water. 
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NOMENCLATURE 


jBj  loss  of  Bf  from  a  part  V,  of  a  soil  column 
b  positive  constant 

C  concentration  of  Bf  in  unfrozen  water 

C+  C/Co 

Ca(0  average  concentration  of  Bf  in  unfrozen  wa¬ 
ter  removed  from  a  part  Vj  of  a  soil  column 
during  the  experiment 

Ct  C^(i)ICo 

Cq  concentration  of  Bf  in  distilled  water  that 

was  mixed  with  dry  soil  to  make  soil  columns 

Cq  concentration  of  Bf  in  unfrozen  water  divid¬ 
ed  by  Co  at  the  beginning  of  the  experiments 

Q  concentration  of  Bf  in  unfrozen  water  for 
segment  i 

C(  solute  concentration  in  the  liquid  phase 

Cs  solute  concentration  in  the  solid  phase 

C  Bf  content,  weight  of  Bf  per  unit  weight  of 
dry  soil 
C+  C/Co 

Cq  initial  Bf  content 

FbCO  mass  flux  of  Bf  from  segment  i  to  segment 

/  +  1 

Fy,(i)  mass  flux  of  water  from  segment  i  to  segment 
/  +  1 


k  effective  segregation  coefficient 

ko  value  of  k  when  u=0 

T  temperature 

Tg  temperature  at  the  cold  end  of  a  soil  column 

Tw  temperature  at  the  warm  end  of  a  soil  col¬ 

umn 

t  time  (days) 

to  22  days 

u  growth  rate  of  a  solid  phase 

Vj  part  of  a  soil  column  consisting  of  segments 
1  through  i 

V  volume  of  each  segment 

Wj  loss  of  water  from  a  part  Vj  of  a  soil  column 

w  content  of  water  (including  ice)  based  on  the 
dry  weight  of  soil 

w*  w/wq 

wq  initial  water  content 

Wi  content  of  water  (including  ice)  of  segment  i 

based  on  the  dry  weight  of  soil 

p  dry  density  of  soil 

pa  average  p  over  all  segments 

Pi  dry  density  of  soil  in  segment  i 

a  standard  deviation 
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Transport  of  Tracer  Br  in  Frozen  Morin  Clay 
in  Response  to  Temperature  Gradients 

YOSHISUKE  NAKANO 


INTRODUCTION 

The  transport  of  solutes  in  frozen  or  freezing  soils 
is  important  for  both  engineering  and  scientific  rea¬ 
sons.  The  presence  of  salts  is  known  to  drastically  al¬ 
ter  the  mechanical  properties  of  frozen  soils  (Sego  et 
al.  1982,  Hivon  1991).  It  is  also  known  that  the  pres¬ 
ence  of  salts  tends  to  reduce  the  rate  of  frost  heave  in 
fine-grained  soils  (Konrad  1990).  The  migration  of 
fertilizer  and  pesticides  during  annual  freezing  is  a 
subject  of  considerable  interest  in  agronomy  (Gray 
and  Granger  1986). 

The  phenomenon  of  freezing  soils  that  contain  sol¬ 
utes  is  complex  and  poorly  understood.  Inasmuch  as 
solutes  tend  to  be  excluded  from  growing  ice,  they 
are  confined  to  the  domains  of  unfrozen  (interfacial) 
water.  The  transport  process  of  solutes  depends  main¬ 
ly  on  the  transport  of  unfrozen  water,  but  it  also  is  in¬ 
fluenced  by  the  diffusion  or  dispersion  of  solutes  and 
complex  interactions  of  solutes  with  the  surfaces  of 
soil  particles  and  ice.  The  amount  of  unfrozen  water 
in  frozen  soils  generally  increases  with  the  increasing 
concentration  of  solutes  (Banin  and  Anderson  1974). 
Since  the  mobility  of  unfrozen  water  tends  to  increase 
with  the  increasing  unfrozen  water  content  (Nakano 
1991),  the  two  transport  processes  of  solutes  and  un¬ 
frozen  water  mutually  depend  on  each  other.  Because 
of  this  mutual  dependence,  the  transport  of  solutes  in 
frozen  soils  is  more  complex  than  that  in  unfrozen 
soils. 

The  rate  of  solute  transport  is  anticipated  to  be 
great  where  the  rate  of  water  transport  is  great.  Since 
the  transport  of  water  is  closely  coupled  with  the 
phase  change  of  water  in  freezing  soils  and  solutes 
tend  to  be  excluded  from  growing  ice,  the  solute  re¬ 
distribution  is  expected  to  be  most  rigorous  near  a 
freezing  front.  Therefore,  the  solute  redistribution 


near  a  freezing  front  has  attracted  the  attention  of  re¬ 
searchers.  Mahar  et  al.  (1983)  investigated  salt  redis¬ 
tribution  during  freezing  of  gravel  saturated  with  a 
salt  solution  of  a  known  concentration.  They  found 
that  a  distinct  freezing  interface  did  not  exist.  Instead, 
a  zone  of  partially  frozen  gravel  with  a  high  unfrozen 
water  content  beginning  at  the  depressed  freezing 
point  isotherm  (freezing  front)  was  observed.  It  was 
also  found  that  the  salt  concentration  in  unfrozen  wa¬ 
ter  increases  toward  the  lower  temperatures  and  that 
the  content  of  salt  in  the  frozen  part  of  gravel  was 
less  than  the  initial  content  (Mahar  et  al.  1983). 

When  the  solidification  occurs  from  the  melt  con¬ 
taining  impurities,  the  composition  of  the  solid  phase 
generally  differs  from  that  of  the  coexisting  melt. 
This  phenomenon,  called  segregation  or  redistribu¬ 
tion  among  researchers  of  crystal  growth,  has  been 
studied  extensively.  Analyzing  the  process  of  redistri¬ 
bution,  Burton  et  al.  (1953)  obtained  a  quantitative 
description  (BPS  theory)  of  an  effective  segregation 
coefficient  k  given  as 

^  =  q/Q=*o/K  +  (i-*oK'’“]  (1) 

where  Cj  and  Q  =  solute  concentrations  in  the  solid 
and  liquid  phases,  respectively, 
adjacent  to  the  solid/liquid  inter¬ 
face 

kQ  =  value  of  k  when  the  growth  rate  u 
of  the  solid  phase  is  zero 
b  =  positive  number  that  depends  on  a 
given  system. 

Weeks  and  Lofgren  (1966)  showed  the  validity  of 
BPS  theory  applied  to  the  problem  of  freezing  salt  so- 


lutions.  By  freezing  columns  of  sand  saturated  with  a 
NaCl  solution  at  constant  rates,  Baker  and  Osterkamp 
(1989)  found  that  BPS  theory  approximately  holds 
true.  Konrad  and  McCammon  (1990)  conducted  step 
freezing  tests  on  a  saturated  clayey  silt  with  various 
salinities.  They  found  that  the  rate  of  cooling  is  the 
main  factor  controlling  solute  redistributions  and  that 
there  is  a  threshold  rate  of  cooling  of  3°C/day  above 
which  no  redistribution  takes  place  (k  =  1). 

The  flow  of  aqueous  solutions  through  unfrozen 
porous  media  has  been  investigated  extensively 
(Nielsen  et  al.  1972).  In  saturated  systems,  the  flow 
of  a  solute  is  described  as  the  sum  of  a  convective 
term  and  a  diffusive  term.  The  diffusive  term  is 
equivalent  to  Pick’s  law  of  diffusion  where  the  diffu- 
sivity  coefficient  includes  the  effect  of  both  molecu¬ 
lar  diffusion  and  hydrodynamic  dispersion.  Such  a 
description  is  not  accurate  for  fine-grained  soils  with 
small  water  contents  because  the  interaction  between 
solutes  and  surfaces  of  soil  particles  becomes  signifi¬ 
cant.  It  has  been  shown  experimentally  that  a  soluble 
salt  and  water  in  which  it  is  dissolved  move  at  differ¬ 
ent  rates  through  fine-grained  soil  (Biggar  and  Niels¬ 
en  1962,  Krupp  et  al.  1972).  For  example,  when  a 
CaCla  solution  containing  ^^ci  and  is  passed 
through  soils,  the  ^®C1  appears  earlier  in  the  effluent 
than  the  ^H.  This  behavior  can  be  explained  by  anion 
repulsion  or  negative  adsorption  of  salt  by  the  nega¬ 
tively  charged  soil  particle  surfaces. 

While  it  has  been  anticipated  that  the  phenomenon 
of  anion  repulsion  is  important  in  the  transport  pro¬ 
cess  of  solutes  in  frozen  soils,  no  supporting  data  are 
available.  The  objective  of  this  study  is  to  examine 
experimentally  the  transport  of  the  Br  ion  in  unsatur¬ 
ated  and  partially  frozen  soil.  The  content  of  Br  was 
intentionally  kept  low  so  that  the  transport  of  unfro¬ 
zen  water  would  not  be  altered  significantly  by  Br . 
This  study  also  aims  to  find  the  effect  of  interactions 
between  Br  and  surfaces  of  soil  particles. 


EXPERIMENTS 

Experiments  were  conducted  in  which  closed  soil 
columns,  initially  with  the  uniform  dry  density  and 
the  uniform  contents  of  water  and  Bf ,  were  subjected 
to  a  constant  and  linear  temperature  gradient.  The  ap¬ 
paratus  used  in  the  experiments  is  shown  in  Figure  1. 

The  apparatus  consists  of  four  parts  made  of  alu¬ 
minum:  an  upper  plate,  a  lower  plate  and  two  end 
plates.  The  upper  and  lower  plates  are  approximately 
23  cm  long  x  20  cm  wide  x  5  cm  thick.  Four  grooves 
with  a  rectangular  cross  section,  20.32  cm  long  x 
1.27  cm  wide  x  0.625  cm  deep,  were  machined  on 
the  upper  surface  of  the  lower  plate.  The  center  lines 
of  these  grooves,  which  are  parallel  to  the  longer  side 
of  the  plate,  are  spaced  3.8  cm  apart.  Grooves  for  O- 
ring  seals  were  also  machined  around  each  rectangu¬ 
lar  groove  in  the  upper  plate  so  that  when  the  upper 
and  lower  plates  are  bolted  together,  four  closed  emp¬ 
ty  spaces  of  rectangular  column  shape  are  formed  to 
be  used  as  containers  of  soil. 

Each  of  the  two  end  plates  has  a  2.54-cm  hole 
drilled  into  it,  through  which  an  antifreeze  mixture 
from  a  temperature-controlled  bath  is  circulated.  The 
two  end  plates  are  positioned  by  four  bolts.  Two 
sides  of  the  upper  and  lower  plates  bolted  together 
have  a  small  taper  that  matches  the  taper  of  the  end 
plates  so  that  the  bolted  upper  and  lower  plates 
smoothly  slide  into  the  spaces  created  by  the  two  end 
plates,  with  assured  close  contact.  The  four  aluminum 
parts,  when  assembled,  are  thermally  insulated  by 
foam  plastic.  A  row  of  10  copper-constantan  thermo¬ 
couples  is  placed  with  equal  spacing  through  the  up¬ 
per  plate  along  the  centerline  of  each  groove  to  mea¬ 
sure  and  record  temperatures  in  a  soil  sample. 

The  soil  selected  for  the  experiment,  called  Morin 
clay,  is  a  marine-deposited  clay  obtained  from  the 
Morin  brickyard.  Auburn,  Maine.  Morin  clay  is  a 
non-swelling  clay  with  a  specific  surface  area  of  60 


2 


m^/g  and  a  specific  gravity  of  2.75  g/cm^.  This  clay 
was  used  in  the  experiments  reported  by  Nakano  and 
Tice  (1990).  For  a  given  dry  density  p  and  initial  wa¬ 
ter  content  wq,  the  exact  amounts  of  oven-dry  Morin 
clay  and  distilled  water  containing  400  pg/g  of  Br  as 
KBr  were  weighed  to  be  packed  into  a  groove.  Mix¬ 
ing  thoroughly  the  weighed  soil  and  water,  we 
allowed  the  mixture  to  set  for  a  few  days  to  attain 
moisture  equilibration.  After  all  four  grooves  were 
packed  with  soil  samples,  we  bolted  together  the  up¬ 
per  and  lower  plates  to  seal  each  of  four  soil  columns 
with  a  rectangular  cross  section  and  froze  them  down 
to  about  -15°C. 

An  experiment  began  by  placing  the  bolted  upper 
and  lower  plates  into  the  space  between  the  two  end 
plates.  A  stable  and  linear  temperature  field  was  usu¬ 
ally  established  within  1  hour  by  adjusting  the  tem¬ 
peratures  of  the  two  baths.  After  a  specified  time 
passed,  the  bolted  upper  and  lower  plates  were 
brought  into  a  coldroom  with  its  temperature  being 
set  at  -10°C.  The  two  plates  were  separated  and  each 
sample  was  quickly  sectioned  into  a  total  of  25  seg¬ 
ments  with  equal  lengths.  Each  segment  was  placed 
in  a  glass  weighing  bottle  and  oven-dried.  The  water 
content  and  the  dry  density  of  each  segment  were  de¬ 
termined  gravimetrically.  The  concentration  of  Bf  of 
each  segment  was  determined  using  the  Dionex  Ion 
Chromatography  Model  2010i. 

RESULTS  OF  EXPERIMENTS 

The  test  conditions  of  the  experiments  are  present¬ 
ed  in  Table  1.  For  all  of  these  experiments,  we  aimed 
to  maintain  the  average  dry  density  of  the  soil  col¬ 
umns  at  about  1.45  g/cm^,  although  the  uniform 
packing  of  soil  in  each  column  was  quite  difficult. 
We  examined  the  uniformity  of  the  soil  columns  by 
using  the  measured  dry  density  of  all  segments  after 
each  experiment.  The  mean  Pa  and  the  standard  devi¬ 
ation  a  of  dry  density  p  are  given  in  Table  1.  Al¬ 
though  these  values  are  affected  by  errors  associated 
with  rapidly  sectioning  a  soil  column  into  25  seg¬ 


ments  with  ideally  equal  lengths,  they  serve  as  an  in¬ 
dicator  of  uniformity  in  packing  and  accuracy  in  sec¬ 
tioning. 

In  Table  1  is  the  temperature  at  the  warm  end 
of  a  soil  column  while  Tg  is  the  temperature  at  the 
cold  end  of  a  soil  column.  For  all  the  experiments, 
the  same  linear  temperature  field  was  applied;  the  du¬ 
ration  of  experiments  was  22  days.  In  the  last  column 
of  Table  1  is  given  the  concentration  of  Br  in  dis¬ 
tilled  water  Cq  that  was  mixed  with  dry  soil  to  make 
the  soil  columns.  In  experiments  1-4,  distilled  water 
without  Br  was  mixed  with  dry  soil,  while  400  pg/g 
Bf  was  used  in  experiments  5-8.  Experiments  1-4 
were  conducted  to  find  the  effects  of  Br  on  the  trans¬ 
port  of  unfrozen  water  by  comparison  with  experi¬ 
ments  5-8,  where  Bf  was  used. 

The  results  of  experiments  5-8  are  presented  in 
Figure  2,  where  nondimensional  quantities  and 

C'*’  are  defined  as 

=wIwq  (2a) 

C^  =  CICq  (2b) 

where  w  is  the  measured  water  (including  ice)  content 
at  the  end  of  experiments,  while  wq  is  the  initial  uni¬ 
form  water  content.  The  unit  of  water  content  is  the 
weight  of  water  and  ice  per  unit  weight  of  dry  soil.  C 
is  the  measured  Bf  content.  The  unit  of  Bf  content  is 
the  weight  of  Bf  per  unit  weight  of  dry  soil. 

^  In  Figure  2  the  nondimensional  quantities  W*"  and 
C’*’  are  plotted  versus  the  segment  number  i  that  is 
assigned  1-25  from  the  warm  end  to  the  cold  end. 
The  linear  temperature  field  imposed  is  also  shown  in 
these  figures.  The  0°C  isotherm  is  located  at  the  in¬ 
terface  between  segments  10  and  1 1  so  that  the  tem¬ 
peratures  of  segments  1-10  are  positive  while  those 
of  segments  1 1-25  are  negative. 

It  follows  from  the  definitions  of  eq  2a  and  2b  that 
w'*'  and  C'*'  would  be  one  if  unfrozen  water  and  Bf 
did  not  move  at  all  during  the  experiments.  It  is  clear 
from  Figure  2  that  both  water  and  Bf  moved  from 


Table  1.  Experimental  conditions. 


Temperature  Temperature 


Experiment 

no. 

wo 

(%) 

Pa 

(g/cm) 

a 

T«, 

rc) 

Tc 

ro 

gradient 

racm) 

duration 

{days) 

Co 

iv-g/g) 

1 

22.5 

1.48 

0.13 

1.00 

-1.50 

0.120 

22 

0.0 

2 

25.0 

1.45 

0.10 

1.00 

-1.50 

0.120 

22 

0.0 

3 

27.5 

1.46 

0.07 

1.00 

-1.50 

0.120 

22 

0.0 

4 

30.0 

1.45 

0.12 

1.00 

-1.50 

0.120 

22 

0.0 

5 

22.5 

1.44 

0.11 

1.00 

-1.50 

0.120 

22 

400 

6 

25.0 

1.46 

0.14 

1.00 

-1.50 

0.120 

22 

400 

7 

27.5 

1.47 

0.13 

1.00 

-1.50 

0.120 

22 

400 

8 

30.0 

1.43 

0.08 

1.00 

-1.50 

0.120 

22 

400 
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a.  Experiment  5. 
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Figure  2.  Nondimensional  contents  of  water 
-1-6  ^^+  and  versus  the  segment  number  ifor 

experiments  5-8. 


the  warmer  to  cooler  parts  in  the  soil  columns.  The 
unfrozen  part  (segments  1-10)  of  soil  columns  lost 
both  water  and  Br .  It  is  interesting  that  a  few  seg¬ 
ments  in  the  frozen  part  near  the  0°C  isotherm  also 
lost  both  water  and  Br  and  that  the  minimum  value 
of  appears  in  these  frozen  segments. 

ANALYSIS  OF  DATA 

We  will  consider  the  balance  of  water  and  Bf  in 
the  part  of  soil  columns  Vj  consisting  of  segments  1- 
i.  We  will  assume  that  the  initial  dry  density  of  each 
segment  is  equal  to  the  final  dry  density  of  each  seg¬ 
ment  measured  at  the  end  of  the  experiment.  The  loss 
of  water  Wi  fi'om  the  part  V,  during  the  experiment  is 
calculated  as 

^  =''Xp7(^o-’^j)  (3) 

where  v  is  the  volume  of  each  segment  and  pj  and  wj 
are  the  dry  density  and  the  water  content  of  segment  j 
respectively.  Similarly,  the  loss  of  Br,  Bj,  from  the 
part  Vj  during  the  experiment  is  calculated  as 

«l=''tp,(Co-Cj)  (4) 

J=1 

where  Cj  is  the  Bf  content  of  segment  j. 

The  calculated  values  of  Wi  for  experiments  1-4 
are  plotted  together  with  those  for  experiments  5-8  in 
Figure  3,  where  open  circles  are  the  data  points  from 
experiments  1-4,  while  solid  circles  are  the  data 
points  from  experiments  5-8.  For  instance,  in  Figure 
3a  the  calculated  values  of  Wj  for  experiment  1  are 
plotted  together  with  those  for  experiment  5.  Experi¬ 
ments  1  and  5  are  conducted  under  the  same  condi¬ 
tions,  except  that  Br  is  absent  in  experiment  1.  Sup¬ 
pose  that  the  soil  was  packed  uniformly  in  each  col¬ 
umn  and  that  Br  does  not  affect  the  movement  of  un¬ 
frozen  water.  Then  it  is  anticipated  that  the  calculated 
values  of  Wj  for  experiment  1  should  be  equal  to 
those  for  experiment  5.  Since  the  uniform  packing  of 
soil  is  quite  difficult.  Figure  3  indicates  that  the  level 
of  Bf  content  used  in  experiments  5-8  does  not  sig¬ 
nificantly  affect  the  mobility  of  unfrozen  water. 

The  calculated  values  of  Wj  and  B,  for  experiments 
5-8  are  presented  in  Figure  4.  It  is  easy  to  see  from 
these  figures  that  both  Wj  and  Bj  attain  their  maxima 
and  that  the  segment  number  of  the  maximum  B,  is 
greater  than  that  of  the  maximum  W,.  Let  Fy,(i)  and 
FbCO  be  mass  fluxes  of  water  and  Bf ,  respectively, 
from  segment  i  to  segment  i  +  1 ,  then  W\  and  Bj  are 
given  as 


'o 


Wi  =  jF^ii)dt 

0 

(5) 

'o 

Bi  =  J  F^ii)dt 

(6) 

0 


where  to  is  the  duration  of  the  experiment  (22  days). 

If  we  assume  that  Bf  is  transported  mainly  by  the 
movement  of  unfrozen  water,  the  mass  flux  Fg  of  Bf 
is  given  as 

F^{i)  =  C;F^{i)  (7) 

where  Q  is  the  concentration  of  Bf  in  unfrozen  water 
for  segment  i.  From  eq  6  and  7  we  obtain 

'o 

B,  =  lc,F^m-  (8) 

0 

For  a  special  case  in  which  Q  remains  constant  at  the 
initial  value  Cq,  from  eq  5  and  7  we  obtain 

Bi  =  CoWi .  (9) 

It  follows  from  eq  9  that  both  W\  and  B;  attain  their 
maxima  at  the  same  segment.  This  obviously  did  not 
occur  in  our  experiments. 

Since  Bf  is  excluded  from  growing  ice  and  is  con¬ 
fined  to  unfrozen  water,  the  concentration  of  Bf ,  Q, 
must  be  a  nondecreasing  function  of  i  for  i  >  10  when 
a  linear  temperature  field  is  established  at  the  begin¬ 
ning  of  experiments.  Because  of  the  Bf  exclusion 
from  ice,  C\  increases  with  i  in  segments  where  ice  is 
present;  as  a  result,  the  segment  number  of  the  maxi¬ 
mum  Bj  becomes  greater  than  that  of  the  maximum 
lYi. 

The  Bf  ion  may  not  be  completely  confined  to  un¬ 
frozen  water.  However,  let  us  assume  this  complete 
confinement,  as  it  appears  to  be  a  good  approxima¬ 
tion.  We  will  also  assume  that  the  phase  equilibrium 
of  water  holds  true  in  the  experiments.  Then,  the  un¬ 
frozen  water  content  is  approximated  by  the  equilibri¬ 
um  unfrozen  water  content  determined  by  the  nuclear 
magnetic  resonance  technique  on  a  separate  sample 
of  Morin  clay.  Under  these  assumptions  the  con¬ 
centrations  of  Bf  in  unfrozen  water  at  both  the  begin¬ 
ning  and  the  end  of  the  experiments  were  calculated. 

The  results  of  calculations  for  experiments  5-8  are 
presented  in  Figure  5,  where  nondimensional  quanti¬ 
ties  Cq  and  C*'(/)  are  defined  as  the  concentrations  of 
Bf  in  unfrozen  water  divided  by  Cq  at  the  beginning 
and  the  end  of  the  experiment  respectively.  For  in¬ 
stance,  the  calculated  values  of  Cq  and  C*‘(t)  for  ex- 
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Figure  4.  Calculated  values  ofWi  and  B 
for  experiments  5-8. 


periment  5  are  plotted  vs.  the  segment  number  i  in 
Figure  5a,  where  solid  and  open  circles  are  the  data 
points  of  Cq  and  C'CO  respectively.  The  values  of  Cq 
are  nearly  one  in  segments  1-14  because  ice  is  absent 
while  the  Cq  increases  with  i  in  segments  14-25  as 
the  content  of  ice  increases. 

The  values  of  C^it)  vary  in  the  unfrozen  part  (seg¬ 
ments  1-10),  but  they  are  not  very  far  from  one.  How¬ 
ever,  the  C^(t)  increases  sharply  from  segments  15-23 
in  comparison  with  Cj  and  as  a  result,  the  values  of 
C^it)  become  twice  as  great  as  the  Cq  values  near  the 
cold  end.  There  is  an  interesting  trend  that  the  values 
of  C^(0  become  less  than  Cq  in  segments  12-17.  The 
behaviors  of  Cq  and  C^it)  described  above  can  be 
seen  also  in  Figures  5b-d,  but  become  less  pro¬ 
nounced  as  the  initial  water  content  wq  increases. 

In  Figure  5  the  values  of  C'(t)  fluctuate  about  one 
in  the  unfrozen  part,  and  the  pattern  of  fluctuation  ap¬ 
pears  to  be  more  or  less  random.  However,  in  the  fro¬ 
zen  part  near  the  0°C  isotherm  there  is  a  common  fea¬ 
ture  that  C+(t)  becomes  less  than  Cq  .  It  is  clear  from 
Figure  4  that  both  water  and  Br  moved  from  the 
warmer  to  cooler  parts  in  the  soil  columns.  The  com¬ 
mon  feature  described  above  indicates  that  unfrozen 
water  and  Bf  may  not  move  at  the  same  rate. 

We  will  introduce  a  nondimensional  quantity  C^ 
defined  as 


II 

(10a) 

Q(0  =  Bi/Wi. 

(10b) 

It  follows  from  eq  10b  that  C^(i}  is  the  (time)  average 
concentration  of  Br  in  water  removed  from  the  part  V\ 
of  a  soil  column  during  the  experiment.  It  is  easy  to 
see  that  the  concentration  of  Br  in  water  removed 
from  Vi  must  be  equal  to  Cq  at  the  beginning  and  C(0 
of  segment  i  at  the  end  of  the  experiment.  The  calcu¬ 
lated  values  of  are  plotted  in  Figure  5  in  compari¬ 
son  with  Cq  and  C^(t). 

In  Figure  5,  we  may  consider  that  the  concentration 
of  Br  in  the  unfrozen  part  did  not  change  significantly 
from  the  initial  concentration  Cq  during  the  experi¬ 
ment.  However,  in  the  frozen  part,  except  for  a  few 
segments  near  the  cold  end,  C^  is  greater  than  both 
Cq  and  C'^{t)  in  all  cases.  Since  the  features  of  Cq, 
C^{t)  and  C^  common  to  experiments  5-8  are  the 
most  pronounced  in  experiment  5,  we  will  examine 
experiment  5  below. 

We  will  consider  a  part  V  of  the  soil  column 
consisting  of  segments  11-18  (Fig.  2a  and  5a).  Since 
C'*’  <  1  in  F  (Fig.  2a),  every  section  of  V  lost  Br  dur¬ 
ing  the  experiment.  Since  w"*"  >  1  in  segments  16-18, 
these  three  sections  gained  water  during  the  experi¬ 


ment.  Segments  11-15  lost  water,  but  C^(f)  is  still  less 
than  Cq  .  As  a  result  C^(t)  becomes  less  than  Cq  in 
every  section  of  V.  The  input  of  both  water  and  Br  to 
V  must  have  come  from  the  unfrozen  part  where  the 
concentration  of  Bf  remained  approximately  at  Cq 
during  the  experiment.  Suppose  that  unfrozen  water 
and  Bf  moved  at  the  same  rate.  Then,  C{t)  must  not 
be  less  than  Cq.  However,  it  is  easy  to  see  from  Figure 
5a  that  C^(0  is  less  than  one  in  segments  1 1-16.  This 
is  obviously  contradictory.  Therefore,  we  may  con¬ 
clude  that  Bf  must  have  moved  faster  than  the  unfro¬ 
zen  water  in  segments  1 1-16  where  C'*'(t)  is  less  than 
one. 

The  features  of  Cq,  C^(0  and  in  experiment  5 
described  above  are  common  to  experiments  6-8.  In 
experiment  6  (Fig.  5b)  C{t)  is  less  than  Cq  in  segments 
11-15,  while  in  experiments  7  and  8  (Fig.  5c  and  d), 
C(t)  is  less  than  Cq  in  segments  12  and  13.  In  Figure  5 
C^  is  greater  than  both  Cq  and  C'*'{t)  in  the  frozen 
part  except  for  a  few  segments  near  the  cold  end. 
These  data  strongly  indicate  that  Bf  moved  faster  than 
unfrozen  water  in  the  frozen  part  of  the  soil  columns. 
In  other  words,  the  repulsion  or  negative  adsorption  of 
Bf  by  the  negatively  charged  Morin  clay  must  have 
moved  Bf  faster  than  unfrozen  water.  Since  the  repul¬ 
sion  of  anions  is  caused  by  the  interaction  between  an¬ 
ions  and  surfaces  of  soil  particles,  the  effect  of  anion 
repulsion  on  the  transport  of  anions  is  expected  to  be¬ 
come  less  pronounced  with  the  increasing  unfrozen 
water  content.  This  general  trend  can  be  seen  in  Figure 
5  with  the  increasing  initial  water  content  wq. 

CONCLUSIONS 

The  movement  of  water  and  Bf  was  measured  in 
unsaturated  and  partially  frozen  clay  columns  subject¬ 
ed  to  linear  temperature  fields.  Both  water  and  Bf 
were  found  to  move  from  the  warmer  to  cooler  parts  in 
the  columns.  The  data  were  analyzed  under  assump¬ 
tion  that  Bf  is  completely  confined  to  unfrozen  water 
in  the  frozen  part  of  the  columns. 

The  concentration  of  Bf  in  water  was  found  not  to 
change  significantly  from  the  initial  concentration  in 
the  unfrozen  part  of  the  columns.  However,  the  data 
from  the  frozen  part  of  the  columns  strongly  indicate 
that  Bf  moved  faster  than  unfrozen  water  because  of 
the  repulsion  of  Bf  by  the  negatively  charged  surfaces 
of  clay  particles.  The  effect  of  anion  repulsion  on  the 
transport  of  Bf  in  the  frozen  part  of  the  columns  was 
found  to  become  more  pronounced  with  the  decreas¬ 
ing  initial  water  content.  It  may  be  concluded  that  the 
anion  repulsion  by  clay  surfaces  plays  a  significant 
role  in  the  transport  of  Bf  in  unsaturated  and  partially 
frozen  clay. 
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